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Recent improvements to OPUCEM, the tool for calculation of the contributions of various models 
to oblique parameters, are presented. OPUCEM is used to calculate the available parameter space 
for the four family Standard Model given the current electroweak precision data. It is shown that 
even with the restrictions on Higgs boson and new quark masses presented in the 2011 Autumn 
conferences, there is still enough space to allow a fourth generation with Dirac type neutrinos. For 
Majorana type neutrinos, the allowed parameter space is even larger. The electroweak precision 
data also appear to favor non-zero mixing between light and fourth generations, thus effectively 
reducing the current experimental limits on the masses of the new quarks, which assume that the 
mixing with the third generation is dominant. Additionally, disregarding the lack of a clear Higgs 
signal from the LHC and focusing only an electroweak precision data comptability, calculations with 
OPUCEM show that, the existing electroweak data are compatible with the presence of a 5th and 
also a 6th generation in certain regions of the parameter space. 



I. INTRODUCTION 



The categorization of the electroweak (EW) corrections based on their contribution types dates back to a study of 
photon propagated four-fermion processes [IJ. In the original nomenclature, the corrections to vertices, box diagrams 
and bremsstrahlung diagrams were all considered as direct" whereas the propagator corrections due to vacuum 
polarization effects were all named as " oblique" since these participate to the computations in an indirect manner. 
An extended review of this approach and its application of the methodology to Beyond the Standard Model (BSM) 
theories helped coining the term oblique parameters [2 usually represented by the letters S, T and U. The Standard 
Model is defined by the values S — T — U — Q for a given top quark and Higgs boson mass. In a BSM theory 
with new fermion doublets, the S parameter estimates the size of the additional fermion sector and the T parameter 
measures the isospin symmetry violation in that sector. In such a model, the U parameter is dependent on W boson 
width, thus insensitive to new physics. 

As the one-loop exact calculations for a model are tedious and error prone, an open-source C/C++ library called 
OPUCEM (Oblique Parameters Using C with Error-checking Machinery) was introduced to improve the reliability 
and reproducibility of computations in scientific publications [T. This library provides functions to calculate the 



* Electronic address: |ece.asilar@cern.ch| 

t Electronic address: echavlanSSgmail.com' 

-^Electronic address: oktay.dogangunSScernxh 
^Electronic address: sinan.kefeli@boun.edu.tr 
^Electronic address: Erkcan. Ozcan@cern. 
**Electronic address: meiimet.saiimmusak .gciu.tr| 
ttElectronic address: Gokhan. Unel@cern.cli 



2 



contributions to the oblique parameters from a number of BSM models. A command line program constitutes an 
example on how to use the library functions and a graphical user interface facilitates the library's use by non- 
programmers. Additionally, tools for plotting error ellipses in the S ~ T plane and extensive documentation are 
provided. Among the physics models currently implemented, one can cite Standard Model with four fermion families 
(SM4)[4i and 2 Higgs Doublet Models (2HDM). Formulas for computing the effect of adding new fermion doublets with 
Majorana-type neutrinos and of the mixing between the quark generations are also available. The oblique parameters 
S, T and U are computed by using both exact one-loop calculations and with some well-defined approximations 
for various models, providing an additional machinery for error checking, beyond the various internal self-crosscheck 
mechanisms. 

Various studies of the SM4 mass and mixing parameters performed using OPUCEM were previously discussed in [5J . 
This note describes the recent additions to the library, and to the command line and graphical tools, as present in the 
OPUCEM version 00-00-07. The overall available parameter space in the SM4 model and the compatibilities of BSM 
theories with further additional generations are discussed. Despite the recent negative results from the Large Hadron 
Collider (LHC) in the search for heavy quarks using about Ifb^^ of data |6], it is shown that there is still available 
room in the SM4 parameter space. The improvements to the OPUCEM package since the last publication [5 , can 
be summarized as follows: (i) the addition of the mixing in the quark sector between the new heavy quarks and the 
second generation quarks, namely the function calls to calculate T value in the presence of the parameter in addition 
to the previously available ^34; (ii) machinery to perform parameter scans through the command line interface, in 
order to facilitate the exploration of the parameter space defined by the fermion and Higgs masses and the quark 
mixings; (iii) extension of the graphical user interface to include the fifth and sixth generation fermion masses for both 
Dirac and Majorana type neutrinos. In the current implementation, if the masses of the N^^ generation fermions are 
set to zero, the calculations are automatically reduced to those for a model of — 1 generations. 

II. OPUCEM ON THE SM4 

OPUCEM and Gfitter are the two open-source tools available for the calculation of the oblique parameters for a 
number of models^ [7]. In addition to SM4, Gfitter's theory repository includes other models such as Extra Dimensions 
and Little Higgs Models. However for the SM4 model, Gfitter does not incorporate the effect of the mixing between 
the quarks of the fourth and other generations, and it is limited only to Dirac type neutrinos. While these two tools 
have been created for different purposes, the fact that two recent implementations of the S, T formulas exist provides 
an opportunity to check their reliability. As the original goal of the OPUCEM library has been to provide well-tested 
and error- free results (as tested multiple times before), a comparison between the results of Gfitter and OPUCEM 
has been performed. For such a comparison, the input reference values and electroweak fit results should be the same 
in both packages. Since the Gfitter results are presented with two different sets of input values (with J7 = forced 
and U = free) , the reference values and the definitions of the error ellipses for the two programs are given in Table |l] 
for the readers convenience. The default values in OPUCEM have been obtained from the LEP Electroweak Working 
Group through private communication (summer 2009 results) and from the Tevatron Electroweak Working Group |ilOj . 
Figure [1] shows a scan of the fourth generation parameter space performed with OPUCEM using the Gfitter reference 
values for the U = free case. This figure can directly be compared to the Gfitter fourth generation parameter scan 
results as presented in [3, as it has been prepared with the same scan and axis ranges for the new fermion and 
Higgs boson masses (as indicated on the figure itself), and with only Dirac-type neutrinos and no mixings between 
generations. The green area is the region of the 5, T plane that is accessible by SM4, as the input parameters are 
scanned through. Red and black ellipses are the 68% and 95% CL contours of the experimentally allowed S, T values. 
The plus sign shows the SM reference point where all the oblique parameters vanish, whereas the blue horizontal and 
vertical lines define the center of the error ellipses given in the first row of Table |T] Comparison of this plot with the 
Figure 12 of [7] shows perfect agreement between the results of the two software packages. 

Table I: Reference values to use for comparison between Gfitter and OPUCEM results, p indicates the correlation coefficient 
between the 5* and T values. 



After the first version of tfiis paper appeared on tfie arXiv, the authors were informed of the G4LHC package 8 that includes the 
software for the SM4 S, T calculations that were performed in [9]. 
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Figure 1: A parameter scan of the SM4 model with OPUCEM using Gfitter default values and parameters. The plot has been 
produced in a format that allows direct comparison with 7 . 

For the rest of this section, the Standard Model studies with four generations of fermions (SM4) is discussed so as 
to determine the effect of the EW precision data on S* — T parameter space. Cases with Dirac and Majorana type 
neutrinos, both stable and unstable, are considered. The aim is to perform a parameter space scan to investigate the 
viability of the SM4 model, given the current status of the searches and the electroweak precision data results. The 
initial values of mass parameters originate from the direct searches, usually summarized by the Particle Data Group 
fPDG)|ill) and updated with publications from current direct search experiments. For the final values a common 
mass limit of 1 TeV, not too far from the partial wave unitarity bound, seems to be preferred choice in the literature. 
Unless stated otherwise, the step size in the mass scans of the S — T plane is taken to be 10 GeV which is a good 
compromise between the execution speed and the accuracy of the results. For the remainder of this work, OPUCEM 
defaults shown in the last row of Table [l] are used in 5* — T parameter space plots. 

A. With Dirac type neutrinos 

For a fourth generation with a Dirac type neutrino, 6 parameters are commonly considered: masses of the new 
fermions, the Higgs boson mass and the sine of the mixing angle between the third and fourth generation quarks. The 
mixing between third and fourth generation quarks, ^34 , plays an important role in the determination of the best 
fit to the EW precision data. Left-hand side of Figure |2] shows the impact of this mixing on the allowed parameter 
space when the PDG mass limits on the new fermions mass values are considered for an unstable Dirac neutrino. In 
this scenario, the initial default mass pattern is set as: m„4 = md4 = 250 GeV, uih = 115 GeV, = 100 GeV 
and m„4 = 90 GeV. Seven different mixing angle values are considered to investigate the effect on the S — T plane: 
I sin034| = 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6. In Figure [2] each scan with a different mixing angle is shown as a contour 
with a different color, the largest contour (hence allowed parameter space) in the S — T plane corresponds to no-mixing 
and the smallest contour, shown with dark blue, corresponds to | sin034| = 0.6. As seen from the intersection of these 
colored contours with the S, T ellipses extracted from the EW precision data, the size of the allowed region gets 
smaller with increasing mixing, such that in order to be compatible with the EW precision data at the 2cr level, the 
mixing between the third and fourth generation quarks has to have an upper limit, |sin034| < 0.6, independent of 
the other parameters. In general, small mixing angles are preferred and so as to keep the fourth generation in the la 
error ellipse, the mixing angle should satisfy | sin 6134 1 < 0.4. 

Since the minimum experimentally allowed mass for a stable neutrino differs from the unstable case, the parameter 
space is re-scanned for the stable Dirac case. As seen on the right-hand side of Figure [2j the decrease in the neutrino 
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Figure 2: Impact of the mixing between third and fourth family quarks on the fourth generation 5*, T space for unstable (left) 
and stable (right) Dirac-type neutrinos obtained using the PDG mass limit values. The color-filled regions show the set of points 
in the plane that is reachable by varying all the fourth generation parameters except the mixing angles. Colors correspond to 
I sin634| varying between 0.0 and 0.6. 



mass limit downs to 45 GeV reflects itself as an expansion of the allowed parameter space with respect to the previous 
case. However, the constraints for the mixing angle remain as before. 

In Autumn 2011, the CMS Collaboration reported new limits on the minimum allowed heavy quark mass values 
corresponding to to„4 > 450 GeV and nidi > 495 GeV at 95% CL |15f[T7] . The parameter space after this update is 
shown in Figurejsj The colored contours correspond to three different mixing angle values | sin034| = 0.0, 0.1 and 0.2. 
In this scenario, due to the reduction of the allowed parameter space, the 2a error ellipse corresponds to a maximum 
value of the mixing angle, | sin034| < 0.2, both for stable and unstable neutrinos. 
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Figure 3: Mixing angle dependence of the fourth generation quarks. S, T parameter space for unstable (left) and stable (right) 
Dirac-type neutrinos using the updated mass limit values from the LHC: mu4 > 450 GeV and > 495 GeV. The 3 colored 
regions correspond to | sinS34|= 0.0, 0.1 and 0.2 (yellow, orange and green, respectively). 

The impact of assuming a stable neutrino is better illustrated in Figure |4] where three different mixing angle values 
for stable and unstable neutrinos are investigated. The considered mixing angle values are |sin^?34|= 0.0, 0.1 and 
0.2 where the allowed regions size decreases as the mixing angle increases. In the plot, darker colors represent the 
unstable Dirac neutrino case and light colors represent the stable case. One can see from the figure that the stable 
Dirac type neutrino slightly enlarges the available parameter space in the negative S direction, independently of the 
mixing angle. 

Results from Summer 2011 which references the CMS outcomes have also excluded a SM Higgs boson in the mass 
range 120 < niH < 600 GeV with 95% CL |12l fT^ . The remaining mass windows for the Higgs boson are for a heavy 
one of mass larger than 600 GeV and a light one between 115 and 120 GeV. Figure [5] indicates these states for the Dirac 
stable neutrino case with | sin 6*34! = 0.0. The light (heavy) Higgs boson scenario is shown with light blue (orange). 
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Figure 4: Mixing angle dependence of stable vs unstable neutrinos with updated mass limit values: mu4 > 450 GeV and 
ma4 > 495 GeV. 



One can see that light Higgs boson is favored over a heavy one, but even a heavy Higgs boson is still compatible with 
a fourth generation according to the EW precision data. 




Figure 5: Light Higgs mass (light blue) and heavy Higgs mass (orange) dependence of the fourth family fermions in stable Dirac 
neutrino case with updated mass limit values: niui > 450 GeV and mdi > 495 GeV. 



B. With Majorana type neutrinos 



A fourth generation with a Majorana type neutrino requires an extra mass parameter, namely the mass of the 
heavier neutrino denoted by ttina- Therefore, seven parameters have been considered for the oblique parameter 
calculations: masses of the new fermions including the mass of the heavier neutrino, the Higgs boson mass and the 
sine of the mixing angle between third and fourth generation quarks. Figure [6] contains a parameter scan for the 
unstable neutrino case, considering both light (left) and heavy (right) Higgs boson cases. The initial fermion mass 
pattern is set as: m„4 = 450 GeV, nidi = 495 GeV, mn = 115 GeV, = 100 GeV, mjv4 = 1000 GeV and m^j^ = 80 
GeV, where the quark mass lower limits are mostly from the CMS results |14l [TBI [T7] . The maximum values of the 
masses have been taken as 1 TeV for all fermions except for the heavier neutrino where a maximum of 4 TeV has 
been considered. The light Higgs mass is tested for two values, 115 GeV and 120 GeV, and the heavy Higgs mass is 
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scanned from 600 GeV to 990 GeV in steps of 30GeV. The mixing angle values under consideration are | sin034| = 0.0, 
0.2, 0.4 and 0.6 represented by contours of different colors: yellow, cyan, orange and green respectively. A simple 
comparison of the two plots reveals that the light Higgs boson is still the preferred scenario, just like it was the case 
for the Dirac type neutrinos. In this scenario, although smaller mixing angles are preferred, the highest considered 
value of I sin 6*34 1 — 0.6 can still be tuned to appear in the Ict error ellipse. On the other hand, for the heavy Higgs 
bosons case, which is less favored, the highest value of the mixing angle allowed by the Icr error ellipse is | sin034| = 
0.2. It can also be concluded that for the Majorana type neutrinos, the additional parameter (mjv4) gives more 
flexibility to the choice of parameters. Therefore Majorana neutrino case provides a wider allowed parameter region 
than Dirac neutrino case. Especially by tuning the ratio of the masses of the new neutrinos, it is possible to make 
some previously unaccessible portions of the parameter space, compatible with the EW precision data. As a flnal 
remark, stable Majorana case was also separately considered but not included in this section. As the only difference 
between the two cases is the decrease of the neutrino mass lower limit down to 40 GeV, all conclusions for the unstable 
neutrino scenario remain also valid for the stable neutrino case, with an increase of the allowed parameter space size. 
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Figure 6: Light (left) and heavy (right) Higgs mass dependence for the third and fourth family fermion in unstable Majorana 
case in terms of different mixing angle values. The colored regions correspond to \smO\ = 0.0, 0.2, 0.4 and 0.6 (Different colored 
contours represent the allowed parameter space are yellow, cyan, orange and green respectively). 



C. Fourth Generation mixing with light generations 

In the literature most experimental results are reported assuming that the fourth generation quarks mix dominantly 
with the third generation, i.e. BF{u4 Wb) = 100% and likewise for the ^4. This choice is partially motivated by 
the gains of using 6-tagging in the data analysis. It is clear that a proper interpretation of the direct quark search 
limits requires consideration of the mixings with light generations. Furthermore, as has recently been shown in [ISj, 
taking quark mixings into account is extremely important in analyzing the EW precision data as it can significantly 
change the favored mass space. 

OPUCEM has been used to test the implications of the quarks mixings in the case of fully degenerate fourth 
generation fermions, which is claimed to have been ruled out by the EW data according to the PDG [TT]. As seen in 
Figure [7] while the fully degenerate case without mixings is indeed ruled out with more than 99.99% CL, it is possible 
to reach better than la agreement with EW precision data for certain choices of the mixing parameters (in agreement 
with the findings of [18]) One should note that since heavy Higgs boson masses tend to increase the S parameter 
while decreasing T, the introduction of the mixings, which increase the T parameter alone, can open up a significant 
portion of the parameter space in general. 

Therefore it is worth repeating some of the parameter scans detailed in the earlier sections taking into account 
the mixings with the light generations. The goal chosen for this section however, is to provide some input to the 
interpretation of the LHC quark search limits, by setting the mass parameters to certain example values and exploring 
the favored values of the mixing parameters. To this end, the values of the mixings allowed by the current measure- 
ments of the 3x3 CKM matrix elements have been determined. It is worth noting that a full treatment of 4 x 4 
CKM with EW constraints is important for determining the overall allowed parameter space, as discussed in |18| . 
However, to perform first order estimations, 3x3 CKM matrix measurements from the PDG and imposing 
the unitarity condition of the 4x4 CKM matrix was sufficient to determine the values of Ki4 = |T^i4p, i = u,c,t. 
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Figure 7: S — T plane is drawn using OPUCEM for the fully degenerate case with the masses of all the fourth generation 
fermions set to 400 GeV. While with no mixing such a scenario is essentially ruled out (Ax^ = 22.1, black point), with mixings 
of I sin024| = 0.07, I sin^34] = 0.14, the Ax^ = 2.0 (red point), which is within la error ellipse. 



Uncertainties on Kn, assumed to be Gaussian, can be obtained by simple error propagation. This procedure yields: 
= 0.0001 ± 0.0006, K2i = -0.10 ± 0.07, K^^ = 0.22 ± 0.12. 
Noting that the allowed space for the magnitude of K14 is much smaller than the others, the focus should be on 
the 2-4 and 3-4 mixings. A large number (0(10^)) of Gaussian-distributed random {K24, pairs were generated 

using the center values and uncertainties computed above. For each generated non-negative K24 and K34, S, T and 
A)(^ from the center of the S-T ellipse were computed, for a Higgs mass of 115 GeV and assuming the neutrinos to 
be of Dirac type. Then using probability as weights, a two dimensional histogram was filled with |sin024| and 
|sin034| values. This procedure was repeated for a number of different fourth-generation fermion masses (to4g) to 
obtain plots of the mixing space as shown in Figure |8] 




LEP EWWG Summer'aS, CL 68.27% 

LEP EWWG Summer'09, CL 95.45% 




Figure 8: 2D histograms of |sinS24| and | sin 634 1, obtained by using probability in the S-T space as weights, for three 
different values of the degenerate fourth-generation fermion mass. (From left to right: m4G = 300 GeV, m4G = 400 GeV, 
W-4G = 500 GeV.) The small black stars at the center of the colored contours indicate the most favored | sin 024 1 - | sin6'34| pairs. 
The scale of the color-axis and the total integral of the weights are normalized to 10^ non-negative (-ft'24, Kit) pairs. 



By searching for the peak in these 2D histograms, the most favored (|sin024|j |sin034|) pair is determined as a 
function of ra4Q. One can observe that the favored 624 decreases slightly as 77140 increases, whereas 8^4 shows a much 
more rapid decrease. It is important to note that even at low masses, favored 024 is not negligible compared to ^34 as 
shown in Figure [9] It is also worth mentioning that the favored values are not far from the results of the preliminary 
analysis presented in |19j . despite relying on much more limited experimental data. Since the branching fraction for 
the decays of the fourth generation quarks into final states involving third generation generation is approximately 
given by BF ~ \Vtb'\^ / {\Vtb'\^ + jV'cfc'P), it is possible to extract the favored value of the BF as well. As seen in 
Figure 10 left plot, this value varies between 90% and 60% as a monotonically-decreasing function of 77i4G- 

An extensive analysis of the scenario with non-degenerate fourth generation fermion masses requires the scan of 
a large parameter space. Therefore, in order to gain some insight, a particular benchmark case is considered in 
which a light Higgs boson decays into a pair of fourth generation neutrinos, avoiding the current Higgs search limits 
from the LHC experiments [201 |2T]. The example benchmark point is taken to be: niH = 130 GeV, m^4 = 60 GeV, 
777^4 = 120 GeV and vary mu4 = 771^4 between 300 GeV and 650 GeV. Such a benchmark case is itself plausible without 
mixing, particularly when 777^ — 777^4 mass difference is large. Therefore the favored values of the mixings come out 
lower than the values in the degenerate-generation case. However, the ratio of the squared mixings, and hence the 
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Figure 9: Favored values of | 810 6*241 (left) and | sin 834] (right) as a function of fourth generation degenerate mass. 
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Figure 10: Left: favored value of BF ~ \Vtb'\^/ {\Vtb'\'^ + |Vc6'|^) as a function of the fourth generation degenerate mass. This 
plot has been prepared a number of times by iterating the full procedure outlined in the text. For a couple of the m4G values, 
we extract slightly different BF values from these iterations. As a measure of the size of the systematic effects in our procedure, 
we plot all such values. Right: same plot for the case of a non-degenerate fourth generation scenario with only degenerate 
quark masses (mq4 = mu4 ~ md4,)- 



favored value of BF, shows very similar behavior to the degenerate case (Figure 10 right plot). 

In conclusion, the fully degenerate case with mixings is checked to be allowed by the EW precision data. The 
favored values of the mixing angles 62A and ^34 tend to imply non-unity branching fractions for the decays of the 
fourth generation quarks into the third generation final states. The favored value of the BF decreases from about 
90% down to about 60% with the fourth generation mass increasing from 300 to 650 GeV. Light Higgs scenarios with 
the Higgs boson decaying to 4th generation neutrinos are allowed by the EW precision data. Since mixing increases 
the T parameter, for heavy charged leptons low values of the mixing would be preferred. However mixing can still 
play a role with behavior similar to the fully-degenerate case. 



III. ON THE FIFTH AND SIXTH GENERATIONS 



In this section, OPUCEM has been applied to the Standard Model with 5 and 6 generations (SMS and SM6) in 
order to show that these are not excluded by electroweak precision data. It is clear that the lack of a Higgs signal 
from the LHC so far implies even stronger indirect limits on a simple extension of the Standard Model with 5 or more 
generations, and might have already excluded some parts of the available parameter space. On the other hand, in 
certain region of the parameter space, the addition of a 5th and 6th generation would reduce the effective cross section 
in come channels even below the SMS values [22 , . Therefore, a study of these extra-generation cases is still interesting, 
both as an input to beyond-SM model building with extra generations, and as an exercise in understanding what the 
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EW data might have indicated had it been taken into account appropriately before. Keeping in mind Autumn 2011 
LHC results for Higgs boson searches, two values of Higgs boson masses have been considered, namely niH = 115 GeV 
and rriH — 600 GeV. The masses of the extra fermions in SMS with Dirac neutrinos are presented in Table |TT] for two 
example benchmark points. The corresponding points in the S — T plane are shown in Figure [TT] left side together 
with the SMS point. It can be seen that all three points are inside the Icr error ellipse. For SMS with Majorana type 
neutrinos, two example benchmark points are listed in Table [Till with their locations on the S-T plane shown on the 
right-hand side of Figure [TT| 



Table II: Two SM5 benchmark points with Dirac type neutrinos. The first point with ran = 115 GeV and | sin 634] = 0.02, 
leading to S — 0.17 and T = 0.17 is on the left and the second point with tuh = 600 GeV and | sin 6*34] = 0.02, leading to 
S = 0.22 and T = 0.22 is on the right. 
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Figure 11: SM3 and four SM5 benchmark points in S — T plane with Dirac (left) or Majorana (right) type neutrinos. The la 
and 2a error ellipses represent the 2009 results of the i7 = fit from the LEP EWWG. The cross corresponds to the SMS with 
niH ~ 115 GeV. In the left plot, the square(star) corresponds to the first(second) SM5 point in Table [III right plot, the 



square(star) corresponds to the first(second) SM5 point in Table III 



Table III: Two SM5 benchmark points with Majorana type neutrinos. First point with rriH ~ 115 GeV and | sin ^34] = 0.07, 
leading to S = 0.11 and T = 0.13 is on the left and second point with tuh = 600 GeV and | smO^^l = 0.07, leading to S = 0.22 
and r = 0.22 is on the right. 
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The SM6 with both Dirac and Majorana neutrinos has also been investigated. Table [TV] contains two example sets 
of mass values for the additional generations with Dirac type neutrinos. Corresponding points in S" — T plane are 
shown in Figure [12] left side. It is seen that although the SM6 points lie further away from the center of the ellipse 
than the SMS point, they are still within 2a error ellipse, thus not ruled out by the EW data. Similarly, two example 
sets of mass values for SM6 with Majorana type neutrinos are presented in Table [V] and the corresponding points in 
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S — T plane are shown in Figure [12] right side. One can see that these points are within Icr error eUipse. It is worth 
noting that the SM6 points with niH — 115 GeV is significantly closer to central value than the SMS point. 



Table IV: Two SM6 benchmark points with Dirac type neutrinos. The first point with niH = 115GeV and | sinS34| — 0.00, 
leading to S = 0.25 and T = 0.26 is on the left and the second point with niH = 600 GeV and | sin 6*34] = 0.01, leading to 
S = 0.31 and T = 0.31 is on the right. 
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Figure 12: SM3 and four SM6 benchmark points with Dirac (left) and Majorana(right) type neutrinos. The la and 2a error 
ellipses represent the 2009 results of the t/ = fit from the LEP EWWG. The cross corresponds to the SM3 with ~ 115 
GeV. In the left plot, the square(star) corresponds to first(second) SM6 point in Table IV In the right plot, the square(star) 
corresponds to first (second) SM6 point in Table [V] 



Table V: Two SM6 benchmark points with Majorana type neutrinos. The first point with tuh ~ 115 GeV and | sin 6*34! — 0.01, 
leading to S = 0.09 and T = 0.12 is on the left and the second point with mn — 600 GeV and | sin 6*1 = 0.05, leading to 5 = 0.20 
and T = 0.21 is on the right. 
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IV. CONCLUSIONS 



Using the enhancements to the OPUCEM library and associated tools, it is shown that given the Autumn 2011 
direct search results from the LHC experiments, a sequential fourth family is still consistent with the electroweak 
precision data. The statement is valid for both cases with Dirac and Majorana neutrinos and also for both light and 
heavy Higgs bosons. It should also be noted that the EW data still prefers a light Higgs for both Dirac and Majorana 
type fourth generation. Consideration of the mixing with light generation quarks enlarges the available parameter 
space, making even the fully degenerate mass scenario compatible with the EW precision data. Such a mixing also 
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reduces the current SM4 quark mass limits, as the favored BF to the third generation is found to be less than 100 
percent. Finally, analysis with OPUCEM shows that a fifth or even a sixth sequential generation has been in the 
realm of the possible, with a lower than the three-generation SM in some example benchmark cases. In conclusion, 
it is worth highlighting that the only way to discover or completely rule out models with additional generations relies 
on further data collection at the LHC and its analysis. 
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